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Abstract 

We consider, in turn, three systems being acted upon by a regularly pulsed 
harmonic potential (PHP). These are i) a classical particle, ii) a quantum 
particle, and iii) a directed line. We contrast the mechanics of the first two 
systems by parameterizing their bands of stability and periodicity. Interesting 
differences due to quantum fluctuations are examined in detail. The fluctua- 
tions of the directed line are calculated in the two cases of a binding PHP, and 
an unbinding PHP. In the latter case there is a finite maximum line length 
for a given potential strength. 

PACS numbers: 03.20.+i, 03.65.-w, 05.40.+j, 74.40.+k 
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I. INTRODUCTION 



Within the modelization of many problems in mechanics, the potential is taken to be 
time-independent, and sets the stage for the ensuing particle dynamics. Quite generally, 
however, the potential may have its own (parametric) dynamics: the stage is shifting beneath 
the actors' feet, so to speak. In quasi-static problems this parametric dynamics will be on 
such slow time scales that it may be safely ignored; or else, an adiabatic treatment may 
be used [|TJ. On the other hand, if the time scale is of exceedingly rapid character, it may 
be easier to forget the potential altogether, and model the parametric variation through 
some fast degrees of freedom. The Langevin equation represents a possible outcome of such 
a procedure There may be situations where neither the quasi-static nor the Langevin 
approaches is appropriate. In these cases, one is obliged to face the time-dependence of 
the potential head-on - this leads to real analytic difficulties, especially in the quantum 
mechanical case, where one must abandon the notion of eigenstates, and tackle the time- 
dependent Schrodinger equation directly. 

The analytic challenge of explicitly time-dependent mechanics is great, and exactly solv- 
able cases are valuable for providing basic insights. The harmonic potential is widely studied 
in many areas of physics, due partly to its inherent solvability. It is also one of the most 
ubiquitous potentials in Nature, due to the existence of near equilibrium states. It is for 
these reasons that the classical and quantum mechanics of the explicitly time-dependent 
harmonic potential have been studied for many years 

One of the most extreme limits within this class of problems is that of the pulsed harmonic 
potential (PHP). That is to say, the potential exists for extremely short instants, between 
which there is no potential whatsoever. This may be represented by 




(1) 



where in particular 




(2) 
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and for regular pulsing r n = nr. (An alternative terminology has arisen in the field of 
quantum chaos ||, in which one speaks of the system being 'kicked' by such a potential.) 
Apart from more obvious manifestations of such a potential (in which an experimentalist 
externally pulses a system with some form of trapping potential), one can envisage such a 
situation arising in the frame of reference of a rapidly moving particle as it regularly passes 
through regions within which a static harmonic potential exists. 

To our knowledge there is no complete treatment of this system. In fact, we suspect that 
its apparent simplicity may have persuaded workers to add complicating features. We are 
aware of similar harmonic models in the field of quantum chaos |],[nj > but these are typically 
considered on a toroidal phase space in order to make closer contact to classically chaotic 
systems. In such a case, it is known that only periodic orbits exist in the quantum case, 
whereas the classical system may be chaotic depending on the curvature of the potential. 
We are content to study the system on the full phase space. Also, a great deal of effort has 
been directed towards the problem of a 'kicked harmonic oscillator' |Tl||, by which is meant 
a static harmonic potential periodically pulsed with some spatially periodic potential. This 
is used to model the behaviour of a trapped particle under the action of a laser. 

To give a unifying flavour to the present work, we shall study three manifestations of the 
problem. In section II we consider a classical particle in a PHP and parameterize the region 
of stability, along with its associated periodic and quasi-periodic dynamics. In sections 
III-V we consider a quantum particle in a PHP. Sections III and IV are concerned with a 
wave packet centered at the origin, which has no classical analogue. Two different analytic 
formulations of the problem are presented, each with their advantages in application. We 
parameterize the stability and periodicity of the dynamics (including a bizarre cycle with 
period r) and also examine the limit of r — > which we compare to the static harmonic 
potential. In section V we examine an off-centered wave packet. Contact is made with the 
classical system via Ehrenfest's theorem, and we also show that the expectation value of 
the Hamiltonian splits neatly into two pieces which are, respectively, purely quantum and 
classical in origin. In section VI we study the third manifestation of a PHP; namely a directed 
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line in thermal equilibrium with a set of planar harmonic potentials. We determine the 
asymptotic transverse fluctuations of the line for binding potentials, and the finite maximal 
length of the line for unbinding potentials. We end the paper with section VII in which we 
give a detailed summary of our results, along with some general conclusions. 

II. CLASSICAL PARTICLE IN PHP 

We consider a classical particle of mass m in a PHP, which is equivalent to the periodic 
impulsive force 

oo 

F(x, t) = -Xx S (t ~ nT ) ■ (3) 

n=l 

The particle will suffer a discontinuous change in momentum with a period of r. In the 
intervening intervals, the particle changes its position with a constant velocity. It is sufficient 
to describe the particle's trajectory by 

x n = lim x{nr — e) , p n = lim p{nr — e) . (4) 

These quantities satisfy the difference equations: 

x n+l = x n + (r/m)p n+ i 

Pn+l =p n - Ax n . (5) 

We must also specify the initial conditions for the particle: x(Q) and p(0). The initial values 
for the difference equations may then be given as: p± = p(0), and x± = x(0) + (r/m)p(0). 

It is convenient to rescale the momentum to p n = (r/m)p n , so that the difference equa- 
tions now take the form 

%n+l — %n + Pn+l 

Pn+l = Pn~ £.X n , (6) 

where 
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£ = Ar/m , (7) 

is the dimensionless coupling to the potential. It is a simple matter to eliminate one set of 
the difference functions ({x n } say), to give the single second-order difference equation 

Pn+2 = PPn+l ~ Pn , (8) 

where (5 = 2 — £. We supply the two required initial data p 1 and p = Pi + — Pi)- 
[One may also proceed by recasting Eq.(||) in matrix form and determining the dynamics 
and stability of the system from the associated eigenvalues. We shall use the second-order 
difference equation in order to make closer contact with the analysis of the quantum system 
in the following section.] 

This difference equation may be easily solved by introducing the generating function 



R(z) = £ z n p n , (9) 



n=0 



which may be inverted via the contour integral, 

Pn = ^ / -Sr*(*) > (io) 

2lTl Jc Z n+L 

where C encircles the origin counterclockwise (with a radius chosen small enough so as not 
to enclose any singularities bar the pole at the origin). 

Summing the difference equation for p with a weight of z n yields the following expression 
for the generating function: 

«(*)= Vt + T - m 

This function has two simple poles located at a± and a<i where 

«i,2 = f ±^ 2 -4)^. (12) 

For |/3 1 > 2 the poles lie on the real axis and one of them has a modulus less than unity. 
Referring to Eq. (PD , we see this implies that \p n \ grows unboundedly with increasing n. 
Thus, stable evolution of the particle is only possible for \(3\ < 2 (corresponding to < 



£ < 4), which we now examine in more detail. Defining a parameter via (3 = 2cos0, we 
have ai 5 2 = e ±l ^. We evaluate the contour integral in Eq.([H]) by noting that the contour C 
may be deformed around the singularities away from the origin such that J c = — J sing . In 
other words, the required integral along C is equal to minus the residues from the two poles. 
Evaluating these residues, and performing some algebraic manipulations, we arrive at the 
result 

Pn+i = [pi sin(n + 1)0 - p o sinn0] . (13) 
sm <p 

It is clear from the above expression that some form of cyclic behaviour with period nr will 
occur when the value of £ (and thus 0) is adjusted so that sinn0 = 0. In this case we have 
= Mir/n, where M = 1, • • • , [n/2], and consequently, p n+ k = (—l) M pk- From Eq.(^) it is 
easy to check that such a value of also implies x n+ k = (—l)xi~. 

Let us classify two types of periodic motion: PMI(n) - a motion for which all measurable 
quantities take on the same values with a period of nr\ and PMII(n) - a motion for which 
only the energy of the system has a period of nr. Now, the energy (between pulses) is simply 
given by 

m fp^ 2 



E„ = — (14) 

Thus, PMII(n) occurs for any value of satisfying sinn0 = 0. However, PMI(n) only occurs 
for values of which are an even multiple of 7r/n. So the simplest PMII motion occurs for 
n = 2, which corresponds to £ = 2. However, the simplest PMI motion occurs for n = 3 
with M = 2, and corresponds to £ = 3. We stress that these periodic motions exist once 
the parameter £ is tuned to an appropriate value, regardless of the initial data (x(0),p(0)). 

As a final remark, we note that there exists one special PMI motion which has a period 
of It. This may be seen directly from the first order difference equations Such a motion 
is possible if one tunes £ = 4 and adjusts the initial data such that p\ = 2x\. 
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III. QUANTUM PARTICLE IN PHP: 1 



In this and the next two sections we shall examine the evolution of a gaussian wave- 
packet under the influence of a PHP. We shall highlight the similarities of the mean motion 
to the classical dynamics described in section II, as well as some subtle effects which are of 
a purely quantum origin. In this section we shall examine the evolution of a wave-packet 
centered at the origin, and we shall use direct evaluation of gaussian integrals to arrive at 
our results. In the next section, we study the same problem, but with the aid of Fourier 
decomposition. Both methods yield the same results, but in surprisingly different formats, 
which are individually suited to the calculation of different quantities. The final section of 
the three is concerned with the evolution of an off-centered wave-packet (which has non- 
zero expectation values for position and momentum and may therefore be compared to the 
classical case). 

Our starting point is Schrodinger's equation for the wave function ip(x,t), with the 
potential given as in Eqs.([l|) and (fj): 

h 2 A 00 

m& = -^ty + if £ S(t - nr) V ■ (15) 

As an initial condition we take a centered gaussian 

i/)(x, 0) = (2B r /tv) 1/4 exp[-Bx 2 ] , (16) 

where B is complex, and B R = Re [5] > 0. This choice is made on the grounds of sim- 
plicity, but we can imagine preparing such a state from the lowest eigenstate of a static 
harmonic potential. The evolution of the wave function between pulses is simply free par- 
ticle propagation. Therefore we can describe the dynamics by the set of functions {i/j n (x)}, 
where 

ipn{x) = lim^)(z, nr — e) . (17) 

e— »0 

It is not a trivial matter to determine the change in the wave function due to a pulsed 
potential. We refer the reader to Ref. for a full discussion of this point (in the statistical 
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mechanics context). The result is that the wave function suffers a discontinuity in phase. 
Thus, the probability density of the particle is unchanged in the immediate temporal vicinity 
of the pulse. With regard to the general pulsed potential given in Eq.([L|), the wave function 
immediately after the pulse is given by 



ipn( x ) = hm^x, nr + e) = ip n exp 



IV [XI 



h 



(18) 



which has a very natural form when viewed from the path intregral perspective |I3] . In the 
present case of a PHP we have 

iXx 2 



ipn(x) = ip n exp 



2h 



(19) 



The free particle propagation between pulses may be written as 



ip n+ i(x) = / dx' G{x — x' , t)i/j^(x') 



(20) 



where the Green function has the familiar form fH 



G (M) = ( TTZTU) ex P 



vmx 

2ht 



(21) 



\2iTiht, 

Combining Eqs. (|1^) and ( p0|) yields the iteration rule (which resembles a transfer matrix in 
the statistical mechanics context) for the functions {ip n (x)}] namely, 

i\x' 2 



tpn+i (x) = J dx' G(x — x', r) exp 



2h 



(22) 



It is convenient to rescale x — > y = x/b, where b = (hr/m) 1 ^ 2 and define the dimensionless 
parameters, £ = Xr/m (cf. Eq.(^) in section II), and rj = 2Bb 2 . The iteration rule now takes 
the form 



M y) = M - 1/2 J dy> exp 
and from Eq . (|T6|) , the initial wave function is given by 



\{y - v'f -%v a 



My') 



(23) 



^(y,0) = ( V R /nb 2 ) 1 /'exp[- V y 2 /2] 



(24) 



where r] R = Re [rj\. 



We shall concentrate on calculating two important physical quantities: the probability 
density P(x,t) = \ip(x,t)\ 2 and the expectation value of the Hamiltonian (between pulses). 
The latter is defined at the moment prior to the pulse: 



E„ 



2m x 



(25) 



but is constant for the duration of the interval between two adjacent pulses. 

The iteration rule (|23[) clearly shows that the wave function will have a gaussian form for 
all times, given that its initial form is chosen to be a gaussian. Thus, we write the general 
form for ip n (in the unsealed x coordinate) as 

„2 



i> n (x) = A n exp 



X 



-a 



n 2b 2 



(26) 



where A n and a n are complex numbers. The probability density is given by 



P n {x) = \A n \ 2 exp[-C n x 2 /b 2 ' 



(27) 



where C n = (a n + cr^)/2 = 7r6 2 |y4 n | 4 , the latter equality following from normalization. From 
Eqs.(|2"5j) and fl2"6] ) one also has 



h \a n \ 2 



2r K + <) 



(28) 



It is useful to define 



I n = I dx i> n {x) = A n {2Tib 2 /a n ) 1 l 2 



(29) 



which may be inverted to yield 



a n = 2nb 2 (A n /I n 



(30) 



Therefore, we may determine all the quantities of interest by evaluating A n (the value of 
the wave function at the origin), and I n (the spatial integral of the wave function). 

As a first step in the evaluation of these two quantities, let us explicitly iterate the 
function tf) n back to the initial condition. Using Eqs. fl23D and (El) we have 
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Mv) = (??>& 2 ) 1/4 (2m)- n/2 J dyn-i ■■■Jdy l 



o exp 



where the n x n matrix M^) has diagonal elements Mqo — V — h - 
diagonal elements = i for \l — m\ — 1, and = otherwise. 
Thus, the wave function at the origin is given by 

where q n = det M*"). Also, integrating the above n-fold integral over y, we find 



I n = (27r6 2 ) 1 / 2 (^/7r6 2 ) 1 / 4 r™/ 2 (-iq n + q^J-W 



(31) 

-i(3, and off- 



(32) 



(33) 



So we may describe the entire dynamics from the set of determinants {q n }- Before explicitly 
calculating these functions, we shall first express the physical quantities of interest in terms 
of {q n }. From Eq.(|30D we may combine the above two expressions to give (for n > 1) 

. fq n + iq n -i\ 



Or, 



(34) 



We therefore have the explicit form for the wave function (for n > 1): 



V„(x) = ( V H /7rb 2 )^ i 



fl/_.i.2\l/4 --n/2 -1/2 



Qn eX P 



ix 2 q n + iq n ~\ 



(35) 



k 26 2 V Qn J) 

The wave function at intervening times may be easily found by propagating the above 



form with the Green function pi]) . Straightforward integration yields 



i/>(x, (n - l)r + Or) = (^/vrfe 2 ) 1 / 4 C n ! 2 [6q n - i(l - ^-x]- 1 / 2 



X6XP W 



g n + ig n _i 



7„ - i(l - 0g, 



n-li 



(36) 



where 6 1 e [0, 1]. 

We may determine the probability density just prior to pulsing ( i.e. the function P n (x)) 
either from our knowledge of A n (cf. Eq.([27|)), or by simply taking the modulus squared of 
ip n as given in Eq. (|35"D . The results of these procedures, although strictly identical, are not 
obviously so, since their equivalence requires the following identity to hold: 
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MQnQn-i] = V R • (37) 

The proof of this statement will be given shortly. The explicit form of the probability density 
is given by (for general times) 

P(x, (n - l)r + Or) = {V . ' J ,. , exp -- . (38) 

\9q n - z(l - 9)q n -i\ { b 2 \9q n -i(l-9)q n - 1 \ 2 } 

The important information contained in the probability density is the temporal evolution 
of the width of the gaussian wave-packet. Just prior to pulsing we have P n ~ exp[— x 2 / , y 2 }, 
where the width 7 n is simply given by 

In = \q n \lo , (39) 

where the initial width 7 = bj\jr\ R ~. Finally, for the expectation value of the Hamiltonian 
between pulses (i.e. the mean energy), we may derive from Eqs.(^), fl34|) and (p7|) the result 

En= 4^jR l^ + ^n-ll 2 - (40) 

We now turn to the evaluation of the set of determinants {q n }. From the definition of 
the matrix M*-™^ it is straightforward to derive the following iteration rule 

q n +2 = -i(3q n +i + q n , (41) 

where (3 = 2 — £. This is very similar to the second-order difference equation (§) that we 
derived previously for the rescaled momenta in the classical system. In the present case, the 
initial data are q = 1, and q\ = r\ — i. Before proceeding to solve Eq. (|41|) , let us first prove 
the assertion ([37]). Multiplying through the difference equation (|4T|) by g* +1 yields 



9n+2^+i = -i(3\q n+i \ 2 + q* n+1 q n ■ (42) 

Now the first term on the right- hand- side is purely imaginary, thus the bilinear object 
R e [9n(2n-i] i s independent of n and therefore equal to 

MQnQn-i] = MliQo] = MV ~ i] = V R , (43) 
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as required. 

The solution of the difference equation ([0]) may be found by use of a generating function 
just as in the classical system. Bounded evolution of \q n \ is again limited to the parameter 
regime \/3\ < 2, which corresponds to < £ < 4. It is convenient to define the parameter 
via 

cos0 = /3/2 = l-£/2, (44) 
as in the classical system. The explicit form for q n is found to be: 

(— i) n 

q n = [iOq — i) sinn0 — sm(n — 1)0] . (45) 

sin0 

It is very tempting to make a direct connection between these results, and those for the 
classical system. However, in doing so we would lose some of the subtleties contained 
within. Consider first, the present quantum dynamics has no classical analogue, since the 
expectation values of the position and momentum remain zero for all times. Second, there 
is a difference between the types of periodic motion in the quantum system considered here, 
and the classical system. In the latter we found it useful to classify periodic behaviour into 
two categories: PMI(n) (all physical quantities having period nr) and PMII(n) (the energy 
having period nr). The condition for motion of the first class was = Mn/n, with M even; 
whereas the condition for the second class was — Mn/n with M odd or even. In the 
present case of quantum motion, it is clear from Eq.(f45|) that periodic behaviour is to be 
expected for sinn0 = which corresponds to = Mn/n with M = 1, • • • , [n/2]. However, 
apart from setting the scale of 0, the integer M plays no other role, since it appears in the 
wave function as a constant phase factor e ~ lM7T / 2 , So there is no distinction between PMI 
and PMII motions in the quantum dynamics of a centered gaussian wave function. 

It is of interest to compare the energy portraits for the classical and quantum motions. 
We scan through values of £ G (0, 2) and follow the evolution of the energy (as given by 
Eqs.flT^I) and (flOD respectively) for twenty or so iterations, all of which values are plotted. 
We refer the reader to Figs.l - 5 for the classical (with typical values of x(0) and p(0)) and 
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quantum (with varying initial parameter rf) portraits. Despite the similarity in the values 
of £ at which the periodic orbits occur, we see a distinct difference in the bounding curves 
Exom{C) and i£ max (£) for the minumum and maximum energies. In the classical case, E min 
is essentially zero and E max is a monotonically increasing function of the coupling £. In the 
quantum case, E min maintains a distinct gap from zero, and the difference E max — E min is 
non-monotonic with £. It would be of interest to determine the analytic properties of these 
bounding curves, but this is beyond the scope of the present work. 

We shall make a more physical connection between the quantum and classical motions 
in section V when we consider an off-centered wave function which allows the particle to 
maintain non-zero expectation values of both position and momentum. First however, we 
present an alternative treatment of the centered gaussian wave function, which will allow us 
to study some other properties of this system more easily. 



IV. QUANTUM PARTICLE IN PHP: 2 

Although the calculations of the previous section were reasonably straightforward, they 
still required the evaluation of cumbersome n-fold gaussian integrals. It might be hoped 
that a simpler derivation of the results is possible, since the wave function evolves by free 
wave propagation between pulses, which may be handled more easily in Fourier space. The 
purpose of this section is to derive an iteration rule for the wave function using this Fourier 
method. The surprise is that the rule (or difference equation) turns out to be non-linear and 



of first order, in contrast to the linear second order rule (|4T|) derived above. The two iteration 
rules must yield the same results of course; but, as we shall see, they are non-trivially related. 
To proceed, we restate that given our initial wave function is gaussian, we can expect 



the wave function at all subsequent times to retain a gaussian form. So, we write 



-if) n {x) ~ exp 



x 2 



(46) 



where we have omitted the prefactor. The wave function at the instant after the pulse is 
given by Eq.(^) which we rewrite here as 

13 



ip+(x) ~ exp 



-{<7n + 



ar 
2P 



(47) 



The next stage of the evolution is free wave propagation, which is most easily expressed in 
Fourier space: ip(k,t + r) = ip(k,t) exp[— ib 2 k 2 /2]. Now the Fourier transform of ip^ takes 
the form 

k 2 b 2 



VvT(A;) ~ exp 



and so 



tp n +i(k) ~ exp 



2(a n + if) 



(48) 



+ i 



2 \£7 n + if 

Finally, inverse transforming the above expression we arrive at 



(49) 



4>n+i(x) ~ exp 



X 

262 



(50) 



where 



+ i . 



(51) 



This iteration rule requires only one piece of initial datum. Evolving the initial wave function 
to t = t — 0, we have 



£7i = 77/(1 +irf) . 



(52) 



Given this non-linear iteration rule, one might expect the system to have some kind 
of non-trivial (chaotic) dynamics. However, with the benefit of hindsight, we know from 
section III that in fact the system has either periodic or quasi-periodic behaviour (in the 
stable regime < £ < 4). This simpler behaviour is not apparent from the non-linear rule 
for {£7 n }, but we can prove it is the case by connecting this rule to the linear second order 
rule for {q n } given in Eq. (f4~l|) . 

Referring to (0), we have the relation (for n > 1) 



a n = -1 



. ( Qn + iq n -i 



Qn 



(53) 
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Substituting this result into the non-linear iteration rule ([51]), we find equality only if q n 
satisfies the linear iteration rule (|4"ID , as required for consistency. In the absence of hindsight, 
the relation ([53]) would be seen as a rather remarkable linearizing transformation. In order 
to probe the general (quasi)-periodic behaviour of this system, the linear iteration rule is the 
description of choice. However, there are two other aspects of this system which are much 
more easily described by the non-linear rule given above. 

The first of these is the existence of a special periodic motion with period r. We may 
identify this by simply demanding that a n +i = &n (= cr) in Eq.([H]). This yields a quadratic 
equation for a with (normalizable) solution 

^=^(4-01 1/2 -|- (54) 

Now, the initial value of a is set by Eq.(^). So this special '1-cycle' is only possible (for 
a given value of £ G [0,4]) for a special value fj of the parameter r\ (which, we remind the 
reader, describes the initial gaussian wave function). We may invert ( p2|) for fj to find 

fj = sin0 + z(l — cos0) = 2 sin(0/2) exp(z0/2) , (55) 

where we have used in place of £, as defined in Eq.([0]). An example of this 1-cycle is 
shown in Fig. 5. This 1-cycle is a purely quantum effect, since the classical system can only 
have cycles of period > 2r for obvious reasons. One might term this motion 'the sound of 
one hand clapping.' 

The second aspect of the quantum system which may be extracted more easily from the 
non-linear iteration is the behaviour of the system as the time between pulses is taken to 
zero. Before taking this limit, it is important to scale out r in other quantities. Now, we 
shall compare the system (in the limit of r — ► 0), to the quantum mechanics of the static 
harmonic potential. For the latter system, the potential is taken to be V s (x) = kx 2 /2. The 
time-averaged potential of the PHP is V(x) = Xx 2 /2t. Thus we shall take the limit r — > 
and A — > with the ratio A/r — > k. It is also convenient to define v n = o n jT. Then the 
wave function just prior to pulsing is given by (cf. Eq. (|46|) ) 
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V>n(;r) ~ exp -— v 



(56) 



From (|5"ip , the iteration rule for t> n takes the form 



1 1 



(57) 



v n + iX/m 



Taking the r — > limit as prescribed above, the above iteration rule tends to the first-order 
differential equation for v (i): 



Now, the ground state wave function of the static harmonic potential is a simple gaussian 



(p(x) ~ exp(— muj s x 2 /2h), where the oscillator frequency u;<( = y/c/m. Referring to the 

differential equation for u(t), we see that t> = u; s is a fixed point of the dynamics. In other 
words, if we initialize the wave function to be the ground state wave function of the static 
harmonic potential, then the wave function will be completely unaffected by the PHP in the 
limit of r — > 0. We stress that there will be no dynamical evolution whatsoever. This is in 
contrast to the wave function of a truly static harmonic potential, which if prepared in the 
ground state, will still have a dynamically evolving phase factor e - Ja -W 2 ^he pulsing of the 
potential essentially resets the clock of the wave function such that the dynamical phase is 
'stuck' at t = 0. This effect may have important consequences for numerical integration of 
the time-dependent Schrodinger equation. If one places the equation on a discrete temporal 
grid, then one is pulsing the potential. From the present example, we see that a pulsed 
potential (on however fine a grid) cannot mimic a static potential. The dynamical phase 
information is irretrievably lost. Subjecting a system to very high frequency pulses was also 
studied recently in the context of controlling decoherence fl5| . 

One can examine the dynamics of the PHP via the differential equation (EH) in more 
detail. For instance, one can examine the evolution of a gaussian wave packet which is not 
tuned to be (p(x). Let us restrict ourselves to an initial condition for ip which is real, and 
therefore completely described by the width of the wave packet, which we denote by Zq, and 




(58) 
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which is related to an effective frequency ui = h/ml^. A straightforward solution of (|58| 
yields 

x 2 ( 1 + (i/r) tan(uty 



ip(x, t) ~ exp 



211 



1 + ir tan(o;t) 

where r = c<j/u; s . Consequently, the probability density P(x,t) ~ exp(— x 2 /7 2 ) with 

/ l + r 2 tan 2 M) \ 
7(t)=/0 ( 1+tan 2 ^) ) ■ 

The expectation value of the kinetic energy is given by 

ftw ^ l + (l/r) 2 tan 2 (^) \ 
W) "T l+tan 2 M) J ' 



(59) 



(60) 



(61) 



In the limit of r — * there is also an effective potential energy (from averaging over the 
many pulses in a given small time interval). The expectation value of the potential energy 
may be found from / dxijj*(KX 2 /2)ip, with ip(x,t) given by (|59|) above. One finds 



hu) (\ +r 2 tan 2 (u;t)~ 
pot(t) = 4^ I, l + tan 2 (^t) , 



(62) 



The total energy is then found to be 



^ot(t) = ^(r 2 + ^ 



(63) 



which is independent of time as expected. 



V. QUANTUM PARTICLE IN PHP: 3 

In this third and final section concerning quantum dynamics in a PHP, we shall consider 
an initial wave function which is off-centered. Thus the expectation values of the position 
and momentum of the particle will be non-zero, and we may make closer contact between 
the quantum PHP and its classical counterpart. We shall utilize the iteration rules from 
both sections III and IV. 

We begin by deriving the first order iteration rules using the method described in section 
IV. Consider an off-centered gaussian at the moment prior to the n th pulse: 
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ip n ~ exp 



x x 

0"n TT7T + d n — 



(64) 



Following the treatment for the centered gaussian (cf Eqs.(46) - (|50|)) we find that the 
iteration rule for {cr n } is unchanged from the form given in Eq.fl5"TD, and the rule for {d n } is 



d-n+i d ri 



0~n+l 



fan + *0 

With the aid of flBlf ) and fl53l) we may rewrite this as 

^n+iQn+i = —id n q n , 
which may be iterated immediately to give the solution 



(65) 



(66) 



d 



n+l 



9n 



(67) 



where do is the parameter introduced to describe the initial off-centered wave function: 
ip(x, 0) ~ exp[— t]x 2 /2b 2 + d x/b]. 

Our present interest will not so much be in {<J n } and {d n }, but rather in the expectation 
values of the position and momentum. Explicitly evaluating the appropriate expectation 
values using Eq.(|64]) we find 

{dn + d* n ) 



x n = dx ip*xip = b- 



(68) 



and 



p n = dx 4>* (-ihd x )t/j 



ih (a*d n - a n d* n ) 



(69) 



b (a n + a*) 

It is convenient to re-express the above relations in terms of {q n } and {d n }; which, due 
to Eq. flBTD , amounts to expressing the expectation values purely in terms of {q n } and do- 
Using the known relations (|37D , (|53|) and (0), we have 

i n b 



X ■■{■ 



2r] J 



:[q n d* Q + (-irq* n d ] 



(70) 



and 
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P» = T^ntilndl + (-l) n q* n d + iqn-id* - z(-l) n <£-Wo] • (71) 



2fop 

We are now in a position to utilize the linear iteration rule (|4lD for the {<2Vi}. Taking 
neighbouring differences between the expectation values of the position and momenta, we 
may use fl4"I|) and the definition 6 = (hr/m) 1 ^ 2 to derive 

7" 

~~ %n — — Pn+1 , (72) 

m 

and 

-p„ = -Ax n . (73) 

A direct comparison between these iteration rules with those given in Eq. @ for the classical 
system reveals that the mean motion of the off-centered gaussian wave function in the PHP 
is identical to the purely classical motion - a non-trivial example of Ehrenfest's Theorem 
TJ. 

Finally, we shall consider the expectation value of the energy. Using the form for the 



wave function given in Eq.(^) we have (cf Eq. 



En = 

2r 



\a n \ 2 (a* n d n - a n d* n ) 2 



K + O K + ® n 



*\2 



(74) 



Referring to Eq. (|28|) , we see that the first term on the right-hand-side is precisely the 'quan- 
tum' energy for the centered gaussian, studied in section III. Referring to Eq. (pP|) , the second 
term on the right-hand-side is determined as nothing more than the classical energy p 2 l /2m 
studied in section II. Thus, the energy of the off-centered gaussian falls neatly into two 
separate pieces: a 'classical' energy determined purely by the expectation value of the mo- 
mentum, and a 'quantum' energy determined purely by the internal fluctuations of the wave 
packet in the co-moving frame. This ends the study of a quantum particle in a PHP. 



VI. DIRECTED POLYMER IN PHP 



Quantum mechanical Green functions may be rewritten as path integrals. It is within 
this formalism that one may appreciate the close mathematical connection between quantum 
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processes and the statistical mechanics of directed lines. We shall use this connection to 
discuss the physics of a directed line in thermal equilibrium with a PHP, which is physically 
realized by a set of static, planar (or sparse ||12||), harmonic potentials. 



We begin by writing the Feynman path integral jn| for the Green function of a quantum 
mechanical particle in a potential V(x, t) (where we restrict our attention to one dimension 
for simplicity): 

y{i)=x 



G(x,t;x ,t ) 



T>y(s) exp 



ds 



to 



m I dy 
~2 [ds 



V(y(s),s) 



(75) 



y(t )=x 

This is to be compared with the (restricted) partition function for a directed line in thermal 
equilibrium with a static potential V(x, z). By 'directed line' we mean a connected path in 
a two dimensional space (x,z), which may adopt any configuration whatever, so long as it 
is directed along the longitudinal (z) direction. The statistical mechanics of such objects is 
of interest in several fields, including directed polymers |16|], and superconducting flux lines 
||17|| . The partition function for a line with one end pinned at (xq, zq) and the other pinned 
at (x, z) is given by 

y(z)=x 



Z(x, z; x Q , z ) = J Vy(s) exp i ~ J ds ^ (jjj^J + 



V(y(s),s) 



(76) 



y(z )=x 

where T is temperature (with Boltzmann's constant set to unity) and k is the elasticity of 
the line. 

The obvious similarity between these two path integrals can be misleading, since it is 
important to remember that the relation between them exists at a strictly mathematical 
level. Given the analytic solution to one path integral, one may infer the solution to the 
other by an analytic continuation from 'real time' t to 'imaginary time' z. However, the 
physical properties of the two systems are generally quite distinct and little may be inferred 
about the physics of one of the systems, if only qualitative information about the physics of 
the other is available. This will become clear in the present case of a PHP, as we shall soon 
sec. 
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The physical meaning of a directed line in equilibrium with a PHP is as follows. The PHP 
itself consists of potentials which only exist on discrete transverse lines, and are regularly 
spaced along the longitudinal axis. These potentials are harmonic and centered at x — 
0. The directed line (which we shall take to be pinned at (0, 0) and to have a length z) 
equilibrates itself in these potentials, meaning that the free energy is minimized as a result 
of the competition between the wandering of the line (entropy), and its elastic and potential 
energies. One may find an application within the field of superconductivity. In an array 
of flux lines, one may construct an approximate (harmonic) 'caging' potential for a given 
flux line, by averaging over the repulsive line-line interactions which it experiences with 
its neighbours. Furthermore, in the strongly layered cuprates |18j (which form the most 
important class of high temperature superconductors) the supercurrents only exist in well 
separated Cu-0 planes. Thus the flux line (meaning the imaginary line connecting the 
planar centres of magnetic flux) will only experience the caging potential in discrete, but 
regular, transverse planes. 

Given the relationship between the quantum system and the directed line at the level of 
path integrals, there is naturally a partial differential equation for Z corresponding to the 
Schrodinger equation. Defining a 'rigidity' v = T/2k, and absorbing T into the potential, 
we have 



The initial condition implicit in the path integral (|76| ) is Z(x,0) = 5(x), but we may gen- 
eralize this to any desired function. Following our earlier work on the quantum system, we 
shall take a gaussian initial condition 




(77) 



The potential is taken to be a PHP, expressed as 




n=l 



(78) 




(79) 
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which would naturally arise from thermal wandering of a directed line from a 5-function 
initial condition. The normalization of Z deserves mention. Whereas in the quantum 
system, the wave function ip is normalized by requiring that / dx\ip\ 2 = 1, the parti- 
tion function has an arbitrary prefactor, as we only require that / dxV(x,z) = 1, where 
V(x, z) = Z(x, z)/ J dx'Z(x', z) is the probability density of the line. 

We shall not enter into any details concerning the analysis of this system, as our results 
may be easily reconstructed from the methods presented in section III for the quantum 
analogue. The present system is described by two dimensionless parameters: an effective 
coupling 



g = 2udg , 



(80) 



and the ratio 



X = ^d/a 2 = {I /a) 



(81) 



which is the square of the ratio of the transverse thermal wandering scale I = (Aud) 1 ^ 2 and 
the initial transverse scale a. The probability density of the line just prior to the n th impulse 
is defined as 



V n {x) = limV(x, nd — e) , 



(82) 



and has the explicit form 



V n (x) 



(1 ~ gw-l/ gn) 
7t/ 2 



nl/2 



cxp 



Qn Qn—1 \ / X 



(83) 



This gaussian form is completely described by one quantity; namely, the width 7„ of the 
probability density defined via V n ~ exp[— x 2 /^ 2 }. Thus 

X 1/2 



ln = l\ 

\q n - q n -ij 

The determinants {q n } satisfy the second order difference equation 



(84) 



q n +2 — I3q n +i — q n , 



(85) 
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where (3 = g + 2, and the initial data are qo = 1 and q± = x + 1- 

Let us first study the case of /3 > 2, which corresponds to an attractive (or binding) PHP 
with g > 0. In this case it is convenient to define a parameter 6 via 

cosh# = (3/2 = 1 + g . (86) 

Then we find 

1 [sinh(n + 1)6 + (x + 1 - 2 cosh 6») sinh n0] . (87) 



sinh 9 

Substituting this solution into Eq.(^4[) yields the final result for the transverse line scale as 
a function of n = z/d : 

1/2 

[X + 1 — cosh^) tanhn^ + sinh# 



In = I 



(88) 



(cosh# — 1)(2 cosh^ — x) tanhn^ + (x ~ 2(cosh0 — 1)) sinh^ 
For the limit of an infinitely long line, n6 — > oo, the function tanhn^ — > 1, and the above 
result simplifies dramatically to 

7oo = - « , (89) 
with finite n corrections ~ e~ 2nd . The dimensionless parameter ( is given by 

C= 1 + sinh^ -cosh^ = 1 - e' 8 . (90) 

Note that for 

C = X (91) 

the transverse scale is asymptotically equal to the initial scale a, meaning that a specially 
tuned PHP can exactly compensate the transverse wandering for arbitrarily long lines. The 
tuned value of the coupling g needed to satisfy Eq. (|9~T|) is 

Since the coupling g must be positive, we see that such a compensating PHP is only possible 
for x < 1, i.e. for a line whose initial scale a > I. This is clear since whatever the PHP 
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strength, the line is free to wander a transverse scale I between pulses, and thus we can 
never restrict the line to 7^ = a if a < I. 

Regardless of the initial transverse scale a, a very strong PHP is expected to strongly 
compress the line to have a scale 700 ~ I. For g ^> 1, one can perform an asymptotic 
expansion on Eq.(^) to find 



7oo = I 1 + — + 



(93) 



The case of a repulsive (or unbinding) PHP is a little more subtle as the line will become 
unstable (meaning the probability density becomes unnormalizable) beyond a critical length 
depending on the strength of the potential. In fact, it is easy to see that for the line to survive 
just one pulse, we require g > —2. It is therefore convenient to define a parameter 9' via 



COS0' = 1 - \g\/2 . 



(94) 



We find for the determinants 



sin 9 



7 [sin(n + 1)9' + ( X + 1 - 2 cos 9') sin n9'\ 



(95) 



and using Eq. 

In = I 



(X + 1 - cos 9') tan n9' + sin 9' 



1/2 



(96) 



cos0' - 1)(2 cos#' - x) tan n9' + (x - 2(cos#' - 1)) sin 9' 
This expression is not valid for arbitrarily large n. There exists a maximum length n*d for 
the directed line, beyond which it no longer exists as a connected elastic structure. The 
value n* may be found by demanding that q n > for all n < n*. Referring to Eq.(p5f) we 
find that n* = [w*], where 



w 



1 

~9> 



* + tan- ( X + 1 7 OSg/ 
2 \ sin 9' 



(97) 



For \g\ — > 2, 6 1 ' — > 7r/2 and 1 < u>* < 2 as expected. 

The more interesting limit of |g| — ► yields the result 

1 



7T 



l^l 1/2 X 



0(\~g\ 1/2 ) 



(98) 
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The transverse scale of the density may also be studied in the limit of small \g\. Referring 
to Eq.(|96|) we find 7 n ~ l\fn for nd' <C 1, which is pure thermal wandering. As n increases 
further, the potential starts to have an effect, and for n9' ~ 7r/2, we find that the transverse 
scale increases linearly with line length 7„ ~ In, with a prefactor depending in a non-trivial 
way on \- 

VII. CONCLUSIONS 

In this paper we have studied in detail the action of a pulsed harmonic potential on three 
systems: a classical particle, a quantum particle, and a directed line. The first and second 
systems share some properties via their mechanics, whereas the second and third share a 
common mathematical basis via the path integral formalism. The pulsing was taken to be 
regular with a period r (or a longitudinal wavelength d in the case of a directed line). 

The classical particle was studied in section II. It was found to have stable (or bounded) 
dynamics as long as the dimensionless coupling £ (cf. Eq. ([?])) lies in the range < £ < 4. 
Otherwise the motion is unstable. For £ < 0, the particle is accelerated to \x\ = oo, whilst 
for £ > 4, the particle 'ping-pongs' with ever increasing amplitude away from the origin. 
In the stable band, there is either periodic or quasi-periodic motion. We introduced two 
classes of periodic motion: PMI(n), for motion where all physical quantities have period 
nr; and PMH(n), for motion where the energy has period nr. We found that the condition 
for PMI(n) motion was that <p = Mir/n where is a convenient parameter defined as 
cos0 = 1 — £/2, and M is a positive even integer satisfying M < [n/2]. Thus the simplest 
PMI motion occurs for n = 3 and M = 2 (corresponding to £ = 3). PMII(n) motion was 
found to occur for = M'lr/n with M' an integer (even or odd) in the range 1 < M' < [n/2]. 
The simplest PMII motion occurs for n = 2 and M' = 1 (corresponding to £ = 2). These 
periodic motions are independent of the initial position x\ and initial momentum p\. We 
also found one special PMI(2) motion which requires £ = 4 and a specially tuned initial 
condition pi = 2mx\/r. 
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The dynamical properties of a quantum particle in a PHP were studied in sections III-V. 
The first two sections concentrated on a gaussian wave packet centered at the origin; whilst 
section V pertained to the case of an off-centered gaussian wave packet, which has non-zero 
expectation values of position and momentum, and may therefore be compared directly to 
the classical particle studied in section II. 

In section III we studied the centered gaussian wave packet in a PHP using the iteration 
properties of the determinants of n-fold gaussian integrals. This led to second order linear 
iteration rules similar to those found in section II. The stability band for the wave packet is 
< £ < 4 as in the classical case. If £ < 0, the wave packet is stretched more with each pulse 
and the width monotonically diverges with time. If £ > 4, the wave packet is squeezed so 
tightly after each pulse that the velocity of expansion of the width of the packet grows ever 
greater after each pulse (due to quantum uncertainty). Within the stable band, there are 
periodic and quasi-periodic dynamics. In the former, the cycles are periodic for all physical 
properties, thus there is no classification into PMI and PMII as in the classical case. The 
condition for a cycle of period nr is <fi = Mn/n with M an integer (even or odd) in the 
range 1 < M < [n/2], where is defined via cos0 = 1 — £/2 as used in the classical case. 
The energy portraits for the classical and quantum particles are shown in Figs. 1 and 2, 
and show several distinct features, especially with regard to the upper and lower bounding 
curves. 

In section IV we analysed the centered gaussian wave packet using Fourier methods, 
which resulted in a first-order, but nonlinear, iteration rule. We showed its equivalence to 
the second-order linear iteration rule of section III. This first order rule allowed two new 
aspects of the problem to be analysed with ease. The first is the existence of a special cycle 
of period r, which is a purely quantum mechanical effect, as a classical system must have a 
cycle of at least 2r. This 1-cycle exists only for an initial complex gaussian wave function 
with an inverse variance tuned to the harmonic coupling via fj = 2 sin(0/2) exp(i0/2). An 
example of a 1-cycle is shown in Fig. 5. The second aspect is the behaviour of the system 
as the time r between pulses is taken to zero. We found that if the initial wave function 

26 



is chosen to be the ground state of a static harmonic potential (with oscillator frequency 
uj = (A/mr) 1 / 2 ), then this wave function is a fixed point of the PHP dynamics as r — > 0. 
This result shows that a pulsed potential with arbitrarily small period cannot mimic a 
static potential, for even though the ground state wave function is a fixed point of the PHP 
dynamics, there is no phase evolution as would be found for a wave function in the static 
potential. The PHP continually resets the phase clock. This causes concern with regard to 
numerical integration of the time-dependent Schrodinger equation, where one may discretize 
time, in which case one is implicitly modelling a static potential by a pulsed potential. 

In section V we allowed the gaussian wave packet to be off-centre, which allowed there to 
be an evolution of the expectation values x and p of the position and momentum respectively. 
Using our previous iteration rules in tandem, we showed that these expectation values obeyed 
the same difference equations as the classical position and momentum, as studied in section 
II, thus verifying Ehrenfest's theorem in a non-trivial setting. Furthermore, we evaluated 
the energy of the wave packet, and found that it split neatly into two pieces: a 'quantum 
piece' equal to the energy of the centered gaussian wave packet, and a 'classical piece' equal 
to p 2 /2m. 

We moved away from the quantum PHP, and in section VI studied the statistical me- 
chanics of a directed line in a harmonic potential which exists only on discrete transverse 
lines which are regularly spaced in the longitudinal direction. This system is the imaginary 
time analogue of the quantum system, as is clear from the path integral formalism. Using 
the same methods as in section III, we found that this system has two qualitatively different 
regimes, depending on the dimensionless coupling g (which is the analogous quantity to £ 
as used in the mechanical systems). For g > 0, we found that the transverse fluctuations of 
the line saturate rapidly to / / \/C, where I is the transverse thermal wandering scale between 
pulses, and £ = 1 — e~ e , where cosh 6* = 1 + g. For very large attractive coupling, the 
transverse scale saturates at I with 0(1/ g) corrections. For g < the line has a maximum 
length n*d beyond which it is destroyed (as a connected elastic entity) by the repulsive PHP. 
An exact expression was derived for n* (as given in Eq. (|97|)), which has the asymptotic form 

27 



for \g\ ^ 0: n* = [w*], with w* ~ n/lgl 1 / 2 . Prior to the line breaking up, the transverse 
wandering grows diffusively for n^l 1 / 2 <C 1, and linearly with n for ^Ig] 1 / 2 — 7r/2. 

These results show clearly the subtle differences which exist between the mechanics of 
the classical and quantum particles in a PHP; and also the differences which exist between 
the statistics of the quantum and statistical path integrals expressions for the PHP. In the 
former case we have seen that the quantum system has another level of complexity beyond 
the mean (or classical) motion. The periodic motion of the centered gaussian wave function 
has no classical counterpart, and its energy portrait certainly deserves more study. In the 
latter case of the quantum versus statistical path integrals, we have seen how no vestige of 
the periodic behaviour of the quantum system remains in the physics of the directed line. It 
has a much simpler asymptotic (n — > oo) behaviour, since all the interesting quantum effects 
are here damped exponentially in n. This serves as a warning that one can only retrieve 
complete quantum information from an imaginary time path integral, if one has complete 
analytic information. 

We believe that these results may also be of some practical interest. It is well known 
how to trap single particles in specially prepared potentials formed from external magnetic 
fields Hl9|| . Thus it is possible to make the trapping potential time dependent by externally 
varying these fields. It would be of interest to use such external variation to mimic a PHP, 
and to test the generality of the results obtained here; namely, the stability band, and its 
associated periodic and quasi-periodic dynamics. 

As detailed in section VI, one can find applications for the directed line in a PHP in the 
field of superconductivity. Of more recent interest in this field is the role of disorder in layered 
materials |[20|| , and its efficacy in pinning flux lines. Disorder is a notoriously difficult effect 
to describe analytically, and there are essentially no analytically solvable cases of directed 
lines in a quenched disorder potential. We consider the generalization of the regularly pulsed 
PHP to one with random pulsing intervals (along with a localized columnar pin) to be a 
prime candidate for such a solvable system, using the framework developed in this paper. 
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FIGURES 
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FIG. 1. Energy portrait for a classical particle of unit mass in a PHP for < £ < 2, and initial 
conditions x(0) = 1 and p\ = 1/2. 
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FIG. 2. Energy portrait for a quantum particle in a PHP for < £ < 2 with h = r = 1, and 
r? = 1 — i/2. 




FIG. 3. As Fig.2 with 77 = 1. 
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FIG. 4. As Fig.2 with 77 = 1 + i/2. 




FIG. 5. As Fig.2 with 77 = 1 + i. Note the 1-cycle at ^ = 2. 
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